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Q-i, 1. Introduction 

^ . The large N limit of four dimensional non-abelian gauge theories is inter- 

esting from the view point of QCD phenomenology 1 and string theory 2 . 
Lattice QCD is a useful technique for extracting fundamental results in the 
large N limit of QCD. Fermions are naturally quenched in the 't Hooft limit 
of large N QCD and this significantly reduces the computational cost in a 
lattice calculation. In addition, there is a concept of continuum reduction 3 , 
namely, physics does not depend on the size of box l A for I > l c and l c is 
a physical critical size. These two observations have been used to extract 
physical results in the large N limit of QCD using numerical techniques on 
the lattice. 



2. Phases of large N QCD 

Large N QCD on a continuum torus I 4 has several phases 4 depending upon 
the size of the torus as shown in Fig. 1. The continuum action has U 4 (l) 
symmetries associated with the Polyakov loops in the four directions and 
the various phases correspond to the number of directions in which this 
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symmetry is broken. The continuum limit is obtained by going to the top- 
right corner of Fig. 1 and different approaches to this corner will result in 
one of the five continuum phases. 

The Oh-phase present for b < 0.36 for all L is an unphysical phase that 
does not survive the continuum limit. The Oh to 0c transition is associ- 
ated with the single plaquette operator opening up a gap around n in its 
eigenvalue distribution. Gauge fields come in disconnected pieces in all the 
Xc-phases due to the presence of the gap in the single plaquette operator 5 . 

The Oc-phase is the confined phase of large N QCD and the lc-phasc 
is the deconfined phase and l c = l/t c - An immediate consequence of con- 
tinuum reduction is that large N QCD does not feel temperatures below 
t c 6 . Numerical analysis has shown that lattice spacing effects are small in 
the Oc-phase and it is sufficient to work at L ~ 9 and N ~ 30 to extract 
continuum results. Therefore, numerical computations can be performed 
on a serial computer and a cluster of computers can be efficiently used to 
generate statistics in a Monte-Carlo calculation. 
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Fig. 2. Extraction of the chiral condensate in large N lattice QCD at a fixed lattice 
coupling. 



3. Chiral symmetry breaking in finite volume 

Since physics does not depend on the box size in Oc-phase, one should show 
that chiral symmetry is spontaneously broken in finite volume in order to 
properly reproduce physics in this phase . The order of limits are impor- 
tant and one has to take the large N limit before taking the quark mass to 
zero at finite phsyical volume. The low lying spectrum of the massless Dirac 
operator shows evidence for spontaneous chiral breaking since the eigenval- 
ues, iX, scales like z = A£iV7 4 with z obeying a universal distribution. The 
chiral condesate, S, is independent of I. Results of a calculation of the chiral 
condensate on the lattice is shown in Fig. 2. Results from chiral random 
matrix theory 8 were used to extract the chiral condensate at a fixed N, 
L and lattice coupling. Two lowest non-zero eigenvalues in the Q = and 
Q = 1 topological sectors were used to show consistency. The plot shows 
that there is a limit as N — > oo and this limit is independent of L. Results 
obtained at different lattice couplings yield ^{V'V'} ~ (0.65) 3 . 
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4. Pions in large N QCD 
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Fig. 3. The pion mass as a function of quark mass in large N QCD. 



Since chiral symmetry is broken in large N QCD even in finite volume, 
one should be able to observe massless pions in finite volume. This result 
emerges in the following manner 9 . Properties of a single quark in a back- 
ground gauge field (x) cannot depend on a shift of (x) — ► A^ (x) + ^f- 
for arbitrary p^ since the U A {1) symmetries associated with the Polyakov 
loops are not broken in the Oc-phase. But the propagator of a non-singlet 
meson will depend on p^, if one quark sees A^(x) — ► A ll (x) + and the 
other sees A tl (x) — > A tl (x) — |y as their respective gauge fields. This is 
referred to as the quenched momentum prescription 10 for the computation 
of meson propagators in the large N limit of QCD. One can use the results 
for the chiral condensate, £(&), and critical lattice size, L c (b), to plot the 
pion mass as a function of m < ■iptp >. The results fall on a single universal 
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curve as shown in fig. 3 and f^l c = , 1 « 0.269. 
5. Chiral symmetry restoration 




Fig. 4. The gap in the quark spectrum as a function of temperature in large N QCD. 



The Oc to lc phase transition is the confinement-deconfincmcnt phase 
transition since the U(l) symmetry associated with one of the Polyakov 
loops is broken. This transition is first order since there is a latent heat as- 
sociated with the single plaqucttc 11 . The fcrmion determinant does matter 
in the lc-phase and it picks the correct boundary conditions for fcrmions 
in the broken direction, namely, anti-periodic with respect to the Polyakov 
loop 12 . The lowest eigenvalue of the Dirac operator with the correct bound- 
ary conditions can be used to study the gap in the lc-phase and one finds 
that chiral symmetry is restored in the lc-phase 12 . Furthermore, the chiral 
transition is first order as shown in Fig. 4. If one were to super-cool the 
lc-phase into the Oc-phase, a second order transition with a square root 
singularity would be observed at 0.93/Z c . 
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6. Phase transition in the Wilson loop operator 

Non-abclian gauge theories in the confined phase are strongly interacting at 
large distances and weakly interacting at short distances. If this is seen as a 
phase transition in some observable, one could use the universal behavior of 
this transition to connect the low energy physics of QCD to the high energy 
physics of QCD. The Wilson loop operator as a function of its size is the 
most likely candidate to study this transition within the Oc-phase of large 
N QCD 13 . Such a transition exists in two dimensional large N QCD 14 and 
it is claimed that the transition in four dimensional QCD is in the same 
universality class. This transition is referred to as the Durhuus-Olesen phase 
transition in Fig 1. 




e/it 

Fig. 5. The phase transition in the expectation value of the Wilson loop operator 
as a function of its size in large N QCD. The distribution of the eigenvalues of the 
smeared Wilson loop operator are compared to the Durhuus-Olesen distributions. Re- 
sults are shown for loops of sizes l/l c = 0.740,0.660,0.560,0.503. They match with 
k = 4.03,2.30,1.41,1.15 where k is the dimcnsionlcss area in two dimension large N 
QCD. k = 2 is the critical area. 

Wilson loop operators suffer from a perimeter divergence in four dimen- 
sions and one has to eliminate it when denning this operator if one were to 
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study its eigenvalue distribution. One possible way to suppress the perime- 
ter divergence is to use smeared operators on the lattice. Numerical studies 
of the eigenvalue distribution of the smeared Wilson loop operator on the 
lattice shows clear evidence for a phase transition as a function of the size 
of the Wilson loop 13 . The universality class is the same as the one found 
in two dimensional large N QCD as shown in Fig. 5 and the critical loop 
size is roughly 0.6l c . 
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